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^ . If the density matrix is treated as an objective description of individual systems, 

it may become possible to attribute the same objective significance to statistical 
^ I mechanical properties, such as entropy or temperature, as to properties such as 

mass or energy. It is shown that the de Broglie-Bohm interpretation of quantum 
^ ■ theory can be consistently applied to density matrices as a description of individual 

On ■ systems. The resultant trajectories are examined for the case of the delayed choice 

interferometer, for which Bell |Bel8?)] appears to suggest that such an interpretation 
is not possible. Bell's argument is shown to be based upon a different understanding 
of the density matrix to that proposed here. 
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1 INTRODUCTION 

^ ! In quantum theory, it is normally assumed that an individual system is described 

^1 by a pure state |$). The density matrix can then enter into the theory in one of 



two ways: 

1. The system is composed of entangled subsystems. Tracing over the Hilbert 
space of a subsystem gives a density matrix which contains the statistics of 
measurements performed only on the remainder of the system: 

piY) = Trx[mX,Y)) mX,Y)\]. (1) 

2. The system is not prepared as an ensemble of identical states, but occurs in 
an ensemble for which the state appears randomly with probability pf. 

P = J2p,\^,) ($,|. (2) 

i 

In both cases, the properties of the density matrix have significant differences 
from an equivalent classical probability distribution. For entangled states, the en- 
tropy of a subsystem can exceed the entropy of the overall state, which cannot 
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happen classically, while for a statistical ensemble, the density matrix does not 
uniquely identify an underlying distribution of pure states. 

For example, the ensembles defined upon a spin-| system, in (jH)), all produce the 
density matrix p = |J, where / is the identity. All the statistical outcomes of an 
experiment can be calculated directly from the density matrix, so there is no possible 
measurement procedure that can distinguish between these three ensembles. 

It has been suggested that the density matrix should be treated as a description 
of a single quantum system rather than as ensemble of systems (see, for example 
|AA98j ). Let us be explicit as to what is meant by this. In Ensemble 1 above, the 
system is in the state |0) half the time. If an observable \A) {A \ is measured on a 
system, the probability of observing the outcome a is 

p{a\0) = \{0 \A)\\ (4) 

Similarly, for systems in the state |1) the probability is p{a\l). It is only over the 
ensemble of states that the result 

p{a) = pop{a\Q) + pip{a\l) (5) 

is obtained. 

If the individual state were described by the density matrix p, then for every 
system the probability of outcome a occurring would be 

p{a) = TT[p\A) {A\]. (6) 

This will not, in general, be equal to either p{a\0) or p{a\l). So the question is 
whether the potential response to a measurement of an individual system must be 
described by a pure state. To distinguish these cases, we will use the symbol g 
to represent the individual state density matrix, and p to represent a statistical 
ensemble from now on. 

As long as the states |0) (0 | and |1) (1 1 occur randomly with probabilities po and 
Pi, it is impossible to establish, by experimental means, any difference between the 
statistical ensemble constructed out of such pure states, and an ensemble of states 
where the statistical outcomes of measurements upon every individual system is 
given by probabilities p{a) that come from a density matrix. This fact alone provides 
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grounds for arguing that it is unreasonable to require, as a matter of principle, 
that individual systems be described by pure states, rather than density matrices. 
The non-uniqueness of the decomposition of the density matrix adds the further 
complication that, even if we are to assume an ensemble of systems is constructed 
from individual pure states, we are unable to determine, by experimental means, 
which set of pure states is involved^. 

A further motivation is provided by thermodynamics. The derivation of thermo- 
dynamic quantities from the statistical mechanics of an ensemble of large numbers of 
systems leads to significant conceptual problems regarding the status of the second 
law of thermodynamics ( |Pop57| ILH.QOl IHPMZ941 lAlbOOl lUffnT] amongst many oth- 
ers). If the density matrix can describe the state of an individual quantum system, 
then that system can have a non-zero entropy 

S = TT[g\n[g\] (7) 

independantly of its belonging to a particular ensemble. This would be expected to 
have a significant affect upon the discussion of the foundations of thermodynamics 
and in particular the discussion of Maxwell's demon and fluctuation phenomena^. 

It should be noted that this suggestion would be more difficult to maintain for 
classical systems. The ontology of classical mechanics describes systems possess- 
ing definite values for all properties. The classical probability distribution uniquely 
defines the underlying states and their distribution, and by a suitably idealised mea- 
surement the particular state can be discovered, non-destructively, in each individual 
case. 

In the quantum case, even for pure states the possession of a particular property 
can only be described through probability distributions. Treating the density matrix 
as having the same ontological status as the pure state then seems to present no 
additional problems. 

However, in |Bel80j . Bell appears to suggest that this is not possible for the de 
Broglie-Bohm interpretation: 

So in the de Broglie-Bohm theory a fundamental significance is given to 
the wavefunction, and it cannot be transferred to the density matrix. 
This is here illustrated for the one-particle density matrix, but it [is] 
equally so for the world density matrix if a probability distribution over 
world wavefunctions is considered. Of course the density matrix retains 
all its usual practical utility in connection with quantum statistics. 

suggesting that only the statistical mechanics of the kind given in )HH96j is valid 
for this interpretation. 

In this paper we show that it is, in fact, possible to apply the de Broglie-Bohm 
approach directly to the density matrix and to construct consistent trajectory solu- 
tions for this. We then apply this to the delayed choice interferometer (Section 0]) 
considered by Bell and show no special problems arise. 

^There are, of course, situations, such as in communication problems |CN 01) . where we are in 
possession of a priori knowledge of the signal states from which a density matrix is composed. In 
such situations the statistical ensemble remains the correct view to take. 

^see jMarn2l Chapter 10] for a brief discussion of this 
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2 BOHM TRAJECTORIES FOR THE DENSITY 
MATRIX 

The de Broglie-Bohm interpretation [Boh52at IBoh52bl IBH93t IHol93j has tradition- 
ally been applied only to pure states. Density matrices have been treated only in 
the context of statistical ensembles [BH96 . To treat the density matrix as a descrip- 
tion of an individual system, we will apply the formalism developed by Brown and 
Hilev |BHOOj . who use the Bohm approach within a purely algebraic framework. 



2.1 The Algebraic Approach 

In |BHnnj , it is suggested that the Bohm approach can be generalised to the coupled 
algebraic equations ^ 



§ = 'lt>.Hl (8) 
4 = -lk,H],. (9) 

Equation ISl is simply the quantum Liouville equation, which represents the con- 
servation of probability, and reduces to the familiar form of 

^ + V.j=0 (10) 

where j is the probability current 

j = R{xr^^ (11) 

m 

in the case where the system is in a pure state g = \ijj) {tp \ and {x \ip) = R{x)e^^^^\ 
The second equation is the algebraic generalisation of the Quantum Hamilton- 
Jacobi, which reduces to 

dt 2m 2mR ^ ' 

for pure states. 

The operator 5" is a phase operator, and this equation can be taken to represent 
the energy of the quantum system. The application of this to the Aharanov-Bohm, 
Aharanov-Casher and Berry phase effects is demonstrated in [BHOOj. 

|BH00j are concerned with the problem of symplectic symmetry, so their paper 
deals mainly with constructing momentum representations of the Bohm trajectories, 
for pure states, and does not address the issue of when the density matrix is a mixed 



[A,B]_ = AB-BA 
[A,B]^ = AB + BA 
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state. Here we will be concentrating entirely upon the mixed state properties of the 
density matrix, and so we will leave aside the questions of symplectic symmetry 
and the interpretation of Equation IHl Instead we will assume the Bohm trajectories 
are defined using a position 'hidden variable' or 'beable', and will concentrate on 
Equation |H1 

The Brown-Hiley method, for our purposes, can be summarised by the use of 
algebraic probability currents 

Jx = Vp{qH) (13) 
Jp = VxigH) (14) 

for which 

^^+[J^,P]_-[Jp,X]_ = 0. (15) 

To calculate trajectories in the position representation (which Brown and Hiley refer 
to as constructing a 'shadow phase space') from this we must project out the specific 
location x, in the same manner as we project out the wavefunction from the Dirac 

ket 'i/j{x) = {x lijj). 

^d{x^^ + {x I [Jx, P]_ \x) - {x I [Jp, X]_ \x) = 0. (16) 

The second commutator vanishes and the first commutator is equivalent to the 
divergence of a probability current 

Vx-J(x)= {x\[Jx,P]_\x) (17) 

leading to the conservation of probability equation 

dP{x) 



dt 



+ Vx ■ J(x) = 0. 



To see the general solution to this, we will note that the density matrix of a 
system will always have a diagonal basis |(;/)a)(even if this basis is not the energy 
eigenstates), for which 

Q = J2'^^ ('^al • (19) 

a 

Note, the Wa are not interpreted here as statistical weights in an ensemble. There are 
physical properties of the state g, with a similar status to the probability amplitudes 
in a superposition of states. 

We can put each of the basis states into the polar form 

i?,(a;)e^^'^(^) = {x (20) 

so the probability density is just 

P{x)=J2WaRaixf. (21) 

a 

The probability current now takes the more complex form 

J{x)=J2waRa{xfVSa{x). (22) 
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So far we have not left standard quantum theory^. We may do this by now 
constructing trajectory solutions X{t), in the manner of the Bohm approach, by 
integrating along the flow lines of this probability current [61393^ .HolQB^ .BHOOj . 
This leads to 

aX(t) ^ J(X(t)) ^ EaWaRa{X{t)yVSa{X{t)) 

dt P{X{t)) EaWaRa{X{t)y ' ^ ^ 

Notice the important fact that, when the density matrix represents a pure state, 
this reduces to exactly the Bohm interpretation for pure states. 

The most notable feature of Equation ESI is that the constructed particle velocity 
is not the statistical average of the velocities {V{t)), that would have been calculated 
from the interpretation of p = J2a |0a) (0a | as an ensemble: 

{Y{t))=Y.wSSa{X{t)). (24) 

a 

This should not be too surprising however. We are interpreting the entire density 
matrix as providing a pilot wave to guide the individual particle motion. All the 
elements of the density matrix are physically present, for a particle at X{t), and 
each state |0a) contributes a 'degree of activity', given by Ra{xY, to the motion of 
the trajectory, in addition to the weighting Wa- If a particular state has a probability 
amplitude that is very low, in a given location, then even if its weight Wa is large, it 
may make very little contribution when the trajectory passes through that location, 
and vice versa. 

Let us consider this with the simple example of a system which has two states 
1 0a) and 1 0b). The probability equations are 

P{X) = WaRaixf + wMxf (25) 
3{X) = WaRa{xfVSa{x)+WbRb{xyVSa{x). (26) 

Let us suppose that the two states |0a) and |0b) are superorthogonal^. This implies 
0a(X)0b(X) ^ for all X. This must also hold for the probability amplitudes 
Ra{X)Rb{X) ^ 0. If the particle trajectory X(t) is located in an area where Ra{X) 
is non-zero, then now the value of Rb{X) ^ 0. The probability equations become 

P{X) ^ WaRa{Xf (27) 
J(X) ^ WaRa{XfVSa{X) (28) 

and so the particle trajectory 

^ VSa{X{t)) (29) 

ot 

follows the path it would have taken if system was in the pure state |0a). In this 
situation, where there is no overlap between the states, then the Bohm trajectories 

^The probability current is a standard part of quantum theory, as its very existence is necessary 
to ensure the conservation of probabihty. 

^Superorthogonahty is defined in |BH93| as a situation where two wavepackets are completely 
non-overlapping in the configuration space of the beable: '^a{X)'^b{X) « 0. It is a much stronger 
condition than orthogonality: (5'a(X) \^i,{X)) — Q. Superorthogonality of measuring device states 
plays a key role in the Bohmian resolution of the measurement problem. 
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behave in exactly the same manner as if the system had, in fact, been in a statistical 
ensemble. 

Now, if we make the assumption necessary to the Bohm interpretation, that 
the initial co-ordinate of the particle trajectory occurs at position X{0), with a 
probability given by P{X[0)), it is apparent that the trajectories, at time t will be 
distributed at positions X{t) with probability P{X{t)). We have therefore consis- 
tently extended the Bohm approach to treat density matrices (and therefore thermal 
states) as a fundamental property of individual systems, rather than statistical en- 
sembles. As we know that the statistics of the outcomes of experiments can be 
expressed entirely in terms of the density matrix, we also know that the results of 
any measurements in the approach will exactly reproduce all the statistical results 
of standard quantum theory. 



3 CORRELATIONS AND MEASUREMENT 

We will now look at how this extension of the Bohm interpretation affects the 
discussion of correlations and measurements. 

The general state of a quantum system consisting of two subsystems will be a 
joint density matrix gi^2- This joint density matrix must be diagonaliscd, before we 
project onto the configuration space of 6o^/i particle positions, using \xi, X2). We can 
represent this projection by a 6 dimensional vector, x, in the configuration space, 
incorporating the 3 dimensions of xi and the 3 dimensions of X2- The probability 
equations are simply 

P{Xi,X2) = '^WaRa{Xi,X2f (30) 
a 

J{xi,X2) = '^WaRa{xi,X2)'^'V:^Sa{xi,X2). (31) 

a 

The probability current can be divided into two 

J{xi,X2) = Ji(a:;i, 0^2) + J2 (2^1, 3:^2) (32) 



where 



Jl{xi,X2) = J2'^aRa{xi,X2YV^^Sa{xi,X2) (33) 

a 

J2(a;i,a;2) = X! ^2)^Vx2S'a(a;i, Xa). (34) 



The conservation of probabihty is expressed as 

dP{xi,X2) 



dt 



+ Vxi • J{xi, X2) + • J{xi, X2) = 0. (35) 



The particle trajectories must be described by a joint co-ordinate X{t) in the con- 
figuration space of both particles, which evolves according to 

dt ~ p{x{t)y ^"^^^ 
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If we separate this into the trajectories of the two separate particles Xi{t) and X2{t), 
this becomes the coupled equations 

gXi(0 _ Ji(Xi(t),X2(t)) 

dt - P{X,{t),X,{t)) ^^'^ 

d^2it) _ J2{X,{t),X,{t)) 

dt - ^^^^ 

We see, exactly as in the pure state situation, that the evolution of one particle 
trajectory is dependant upon the instantaneous location of the second particle, and 
vice versa. 

The first special case to consider is when the density matrices are uncorrelated 

gi,2 = Qi® 92- (39) 
The probability equations reduce to the form 

X2) = Pixi)P{x2) = WaRa{xif ^1 WbRb{x2f (40) 

a b 

J{xi,X2) = P(x2) Ji(xi) + P(a;i) J2(x2) (41) 

where 

) (42) 

a 

J2(a:2) = J2^bRbix2f'^^2Sbix2). (43) 

b 

The resulting trajectories 

dx^ _ Mx.jt)) 

dt - P{X,{t)) ^^^^ 

dx^ _ UX2{t)) 

dt - P{X2{t)) 

show the behaviour of the two systems are completely independant. 
Now let us consider a correlated density matrix 

1 

Qi,2 = 2 I + \^bXb) {(pbXb I) (46) 

where the |0) states are for system 1 and the \x) states are for system 2. The polar 
decompositions 

Ra{Xl)Ra{x2)(^'''^'^-^'''^'^ = {Xl, X2 {(^aXa) (47) 

i?f,(xi)i?f,(a;2)e'^''(^i)+^''(^^) = (xi,X2|te) (48) 
lead to probability equations 

Pixi,X2) = ]^(Ra{xrfRa{x2f + Rb{xifRb{x2f) (49) 
i{Xu X2) = \ {Ra{XlfRa{x2f{y^^Sa{Xi) + V ^^Sb{x2)) 

+ Rb{xifRb{x2f{V^^Sb{xi) + V^^Sb{x2))) . (50) 



The trajectories, X{t), are then given by 

aXi(t) Ra{Xi{t)fRa{X2{t)yVx,Sa{Xi{t)) + Rt{Xi{t)f R,{X2{t)fV St{Xi{t)) 

Tf) = 

dt Ra{X,{t)yRa{X2{t)y + R,{X^{t)YR,{X2{t)f 

(51) 

dt Ra{x,{t)yRa{X2{t)y + Rb{x,{t)yR,{X2{tw 

(52) 

Now in general this will lead to a complex coupled behaviour. However, if either 
of the states \(f)) or \x) are superorthogonal, then relevant co-ordinate, Xi or X2 
respectively, will be active for only one of the Ra or R^, states. For example, suppose 
the Ix) states are superorthogonal 

Ra{X2)Rb{X2) ^ 0. (53) 

For a given location of X2, only one of these probability densities will be non-zero. If 
we suppose this is the \xa) wavepacket, then Rf,{X2)'^ ~ 0. The trajectory equations 
become 

<9Xi(t) R,{X,{t)fRa{X2{t)yVj,,Sa{X,{t)) 

771 = 

dt Ra{x,it))maiX2{t)y 

= Vx,5a(Xi(t)) (54) 

9X2(t) R^{X,{t)fRa{X2{t)fV^,Sa{X2{t)) 



m- 



dt Ra{x,{t)yRaiX2{t)y 

= Vx.>?a(X2(t)). (55) 

Both trajectories behave as if the system was in the pure state \(f>aXa)- If the 
location of X2 had been within the \xb) wavepacket, then the trajectories would 
behave exactly as if the system were in the pure state \(l)bXb)- The trajectories, as a 
whole, behave as if the system was in a statistical mixture of states, as long as at 
least one of the subsystems has superorthogonal states. 

The loss of phase coherence does not play a fundamental role in the Bohm theory 
of measurement. It is the superorthogonality that is important. This carries directly 
over into the density matrix description. It is a simple matter to generalise the above 
arguments to a general N-body system, or to consider states where the diagonalised 
density matrix involves entangled states. 



4 THE DELAYED CHOICE INTERFEROME- 
TER 

We will now consider the Bohm trajectories for the delayed choice version of the 
two slit experiment. The basic experimental arrangement is shown in Figure 

The Bohm trajectories for pure states in this interferometer have been discussed 
extensi vely ^ESflW92l IDH!t981 IDF(4Z93l IESSW98L WTM ITCTsl LScu98[ IHHMnnl 
IMarn2j . in the context of the claim in jESSW92] that the trajectories are "surreal- 
istic" . 
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Figure 1: Delayed Choice Interferometer: interference (left); "which path" (right) 




Figure 2: "Which Path": measured (left); not measured (right) 



We will not revisit these arguments here but simply note that, when a measure- 
ment of the path is made while the particle is in the interferometer (by detectors 
and Cd) the Bohm trajectories follow the paths like Figure ISfleft). In the absence 
of such a measurement, the Bohm trajectories for the particle are deflected in the 
interference region R (Figure right)) so 'fooling' the delayed choice of 'which path' 
made by detectors at Di and D2. ^ 

For the density matrix we will be considering the same experimental arrange- 
ment, but the atomic state entering the arms of the interferometer (after the beam 
splitter region B) is the mixed state 

^ {\<j)u{x, h)) (0„(x, U) I + |(/)d(x, h)) {(t)d{x, h) I) (56) 

where represents a pure state in the upper arm of the interferometer and 0^ is a 
pure state in the lower arm. No interference effects are expected in the region R. 
We will describe the Bohm trajectories for this in the cases where: 

1. The mixed state is a physically real density matrix q\ 

2. The mixed state is a statistical mixture p; 

3. The mixed state is a physically real density matrix, and a measurement of the 
atomic location is performed while the atom is in the interferometer. 

^For a critique of the 'one bit detectors' critical to the analysis of jESSW92] see jMar02l Chapter 

3] 
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4.1 Physically real density matrix 



While the atom is in the arms of the interferometer, the wavepacket corresponding 
to and that corresponding to \(t)d) {4>d\ are superorthogonaL The trajec- 

tories in the arms of the interferometer are much as we would expect. However, 
when the atomic trajectory enters the region i?, both wavepackets start to overlap. 
The previously passive information in the wavepacket from the other arm of the 
interferometer becomes active again. 

No interference fringes occur in the region i?, and if phase shifters are placed in 
the arms of the interferometer, their settings have no effect upon the trajectories^. 
However, the trajectories do change in R. The symmetry of the arrangement, and 
the 'no-crossing principle' for the flow lines in a probability current, ensures that 
no actual trajectories can cross the center of the region R. The Bohm trajectories 
follow the 'surrealistic' paths similar to those in Figure right), even in the absence 
of phase coherence between the two arms of the interferometer. 



4.2 Statistical Ensemble 

We have seen that, even in the absence of phase coherence, the Bohm trajectories 
for the density matrix show the surrealistic behaviour. Does this represent an un- 
acceptable flaw in the model? To answer this, we now consider the situation where 
the density matrix is a statistical ensemble of pure states. This situation should 
more properly be described, for the point of view of the Bohm approach, as a finite 
assembly^ of many individual systems. 
First consider the assembly 

Pi = |0ai) (0ai| ® l^aa) (^aal ® l^as) ((^asl ••• 

= Iii\K){<PaA (57) 

where Ui = u ot d with a probability of one-half. As the assembly consists entirely 
of product states, the behaviour in each case is independant of the other cases. 

If the state is then the trajectories pass down the u-branch, and go 

through the interference region without deflection. Similarly, systems in the \(pd) {4'd I 
state pass down the d-branch and are undefiected at R. These trajectories are what 
we would expect from an incoherent mixture. 

However, now let us consider the assembly 

P2 = n,|0,j (58) 

where bi = + or — occur with equal probability and 

\M = ^m + m) (59) 
= -^m-m)- (60) 



^To observe interference fringes we would need a density matrix that diagonalises in a basis 
that includes non- isotropic superpositions of \(f)u) and \4>d)- 
®The term "assembly" is taken from ^PerQ.Sj . 
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This forms exactly the same statistical ensemble. Now, however, in each individual 
case there will be interference effects within the region R, it is just that the combi- 
nation of these effects will cancel out over the ensemble. If we were to measure the 
state in the (+, — ) basis, then we would be able to correlate the measurements of 
this to the location of the atom on the screen and exhibit the interference fringes. 
The Bohm trajectories for the assembly p2 all reflect in the region R and display 
the supposed 'surrealistic' behaviour. 

There are no observable consequences of the choice of the different assemblies to 
construct the statistical ensemble^. Consequently, if we are only given the density 
matrix of a statistical ensemble, we are unable to say which assembly it is constructed 
from and cannot simply assume that the underlying Bohm trajectories will follow 
the pattern in Figure [U It is only legitimate to assume the trajectories will pass 
through the interference region undeflected if we know we have an assembly of \(pu) 
and \(pd) states, in which case the Bohm trajectories agree. Thus we conclude the 
behaviour of the trajectories for the physically real density matrix cannot be ruled 
out as unacceptable on these grounds. 

4.3 Measuring the path 

Finally, we consider what happens when we have the physically real density matrix 

^ = ^ h)) {4>u{x, h) I + \<j)d{x, h)) {<j)d{x, ti) I) (61) 

and we include a conventional measuring device in the u-path. The measuring device 
starts in the state \^q). If the atom is in the state \(pu), the measuring device moves 
into the state The states |^o) and are superorthogonal. 
If we now apply the interaction to the initial state 

Q^\^o){^o\ (62) 
the system becomes the correlated density matrix 

{<Pu^i\ + \Mo) {<Pd^o\). (63) 

As we saw above, as the measuring device states are superorthogonal, the system 
behaves exactly as if it were the statistical ensemble. This is true even when the 
atomic states enter the region R. The Bohm trajectories of the atom pass undeflected 
through in the manner of Figure left). 

We conclude that the Bohm trajectories for the density matrix cannot be con- 
sidered any more or less acceptable than the trajectories for the pure states. 

5 CONCLUSION 

By extending the de Broglie-Bohm interpretation to cover density matrices, we 
showed it was possible to consistently treat the density matrix as a property, not of 
an ensemble, but of an individual system. 

^It is interesting to note that if we were to measure the assembly pi in the (+, — ) basis we 
would still obtain interference fringes! 
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It is worth asking why Bell seemed to suggest that this was not possible. An 
examination of |Bel80j shows that the density matrix Bell considered arose, not as 
a fundamental description of an individual system, but through tracing over part of 
an entangled system in an overall pure state, as in Equation 

Bell considers a pure state that interacts with a measuring device, in an arrange- 
ment similar to Figure |2lJ^ left). After the detectors have interacted with the atomic 
beam, the full state is: 

^{Huii) + \ct>dio)). (64) 
Tracing over the detector states leaves the density matrix 

P=\i\4>u) + (0d|). (65) 

Bell wishes to construct trajectory solutions which resemble Figure |21^1eft) based on 
the density matrix for the atomic beam alone. This fails, as the entangled state 
trajectories are non-locally dependant upon the location of beable of the measuring 
device. 

The symmetry of Equation IHSl and the 'no-crossing' principle ensure that a tra- 
jectory solution based upon the reduced density matrix alone must produce tracks 
as in Figure |2f right). It is only the superorthogonality of the detector states 
that allows the atomic beam states to act as if they are a statistical ensemble of 
states and cross over in the interference region^''. 

However, the density matrix we have considered in this paper does not arise 
either as a reduced density matrix, or as a statistical ensemble of pure states. It is 
the complete description of an individual system. For this, we have shown that a 
consistent and coherent trajectory interpretation in the manner of the de Broglie- 
Bohm interpretation is possible. 
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